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Abstract 

We use the Kondo lattice model to investigate the possibility of ferromagnetism and half- 
metallicity in local moment systems. Using the spectral density approach and making use of the 
fact that the spontaneous magnetization of local moment and the itinerant electron polarization 
Ch ' are coupled, we derive an expression for the paramagnetic susceptibility. The magnetic ordering 

o : 

CJ . temperature is determined from the singularities of the susceptibility. The magnetic phase 

diagram is constructed in T — ra(band filling) plane. It is found that ferromagnetism is possible 

only for small values of n. It is also found that the temperature drives the transition of the 

^ . system from half-metal to metal. 
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1 Introduction 

Local moment systems constitute one of the most fascinating fields in both theory and experi- 
ment in condensed matter physics^'^. In general there is a consensus that local moments play a 
decisive role in a variety of phenomenon such as magnetism, heavy fermions, high temperature 
superconductivity, colossal magnetoresistance and spintronics. Many of the real systems that 
are classified as local moment systems consist of both the localized and itinerant electrons. 
In the case of rare earth compounds and alloys, for example, the electrons in the partially 
filled 4f-shell are highly localized and maximize the magnetic moment according to Hund's 
rule and thus constitute the localized moment. The electrons in 5d6s-bands constitute the 
itinerant electrons. Other prominent groups of materials whose electronic structure somewhat 
resembles that of rare-earth materials are the 3d transition metal oxides ^ like Lai^xSr^MnO^ 
or Lai-xCaxMnOs and magnetic semiconductors ^'^ such as Ga{Mn)As. Based on the spin- 
resolved photoemission studies ^ it is identified that the manganese perovskite Lao.TiSro.sMnOa 
is a half-metallic system well below the Curie temperature. The similarity to the rare-earth 
compounds stems from the fact that the 3d-band is split by the crystal field into the doubly 
degenerate eg- and the triply degenerate t2g-bands. The t23-band is occupied by three electrons 
with their spins aligned parallel due to the strong Hund's rule coupling and this provides the 
localized spin back-ground with S = 3/2 in which the itinerant e^-band electrons move. These 
electrons again couple ferromagnetically via Hund's exchange coupling with the localized spins. 
Theoretically, the basic model for understanding the magnetic phenomena in systems where 
local magnetic moments couple ferromagnetically to itinerant carriers is the Kondo lattice 
model(KLM). The KLM in the strong coupling limit is also known as the double exchange 
model. Since the many-body problem of the KLM is not exactly solvable, several approxi- 
mation schemes such as the coherent potential approximation^, the second-order perturbation 
theory^, the moment conserving decoupling approach ^, the dynamic mean-field theory^"^^ and 
the interpolation self-energy ansatz^^ have been developed. Even though all these schemes start 
from different limiting cases of the model, the conclusions drawn from them show remarkable 
similarities. In addition to these approximate methods, spectral density approach (SDA) has 
also proven to be very convenient to study various many-body problems. SDA is a reliable 
approximation in the strong coupling limit so long as the one-electron spectral density has 
predominantly a two-peak structure. In the literature SDA has been successfully used for the 
investigation of both the attractive^*^ and the repulsive^^ Hubbard model, the periodic Anderson 
model^^ and the t — J model^'^. Very recently this approach was also used to investigate^^ the 
interacting spin waves in the KLM. In this work the KLM is mapped on to an effective Heisen- 
berg model and then SDA is applied. In addition, SDA was also applied to repulsive Hubbard 



model extended by Jahn- Teller interaction, to study^^ the interplay between magnetism and 
structural distortion in strongly correlated electron systems. 

In the present work, we attempt an approximate solution of the KLM using the SDA. The 
aim of the present investigation is to study the conditions under which the interband exchange 
and band filling together cause a collective (ferromagnetic) ordering and half-metallicity of the 
local-moment system. 

The present investigation is organized as follows. In the next section we introduce KLM 
and discuss the approximation used to solve the underlying many-body problem using SDA. 
To determine the magnetic ordering temperature, we derive the expression for static magnetic 
susceptibility in the third section. The fourth section is for a discussion and interpretation of 
the results. In the last section we conclude our results. 

2 Model and theory 

Magnetic materials whose magnetic properties depend on a system of localized electrons which 
are coupled via exchange interaction to itinerant conduction electrons are described by the 
KLM, which is given by 

H = Hs + Hsf (1) 

Hg describes the itinerant electrons in the presence of a homogeneous magnetic field B. 

Hs = Y^ (Tij - Z^lXBB5ij) cl^Cja = XI (^k - Zal^BB) 4^Ck^ (2) 

Tij are the hopping integrals being connected by Fourier transformation to the free Bloch 
energies ek. /is is the Bohr magneton. c|^(cjo-) and cj^^(cko-) are the creation(annihilation) 
operators of itinerant electron with spin a {a =t, J,) at lattice site Rj and with a wavevector k. 
Zrj is a sign factor, z^r = S^j^ — 6„i. The second term in Eq.(l) is the interband exchange term 
with coupling strength J, written as intra-atomic interaction between the conduction electron 
spin ai and the localized magnetic moment represented by the spin Si. 

j 
It is more convenient to express this exchange interaction in the following second quantized 
form: 

Hsf = -\jY.i^.S^A^c,„ + Src]_,c,^) (4) 
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Where 5/ = S^ + iz^S^. 



We are mainly interested in the itinerant electron properties. So the relevant quantities to 
be calculated are the single-electron Greens functions and the spectral function. 



GUE) = ((ck.; cL))^ ; S^E) = -^'^GUE + tO' 



(5) 



In the atomic limit and for classical spin, this Greens function has two poles at ±JS/2. There- 
fore it is realistic to make a two pole ansatz for spectral density in the finite band width case 
also: 

Si^^iE) = Y.a^{k,j)6{E - E,{k,j)) (6) 

i=i 

Here the spectral weights ^^(k, j) and the poles Ecr(k,j) are unknown and are to be determined. 

For this, we use two definitions of spectral moments: 
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The two-peak ansatz for the spectral density involves only four unknowns. Therefore it is 
sufficient to calculate only four moments. Calculating the commutators in Eq.(8) and using 
Eq.(6) in Eq.(7) we get the set of equations to determine the spectral weights aa(k,j) and the 
poles E„(k,j) as follows: 

a^(k,l) + a^(k,2) = 1 



E^{k, l)a,(k, 1) + E^{k, 2)a,(k, 2) = ek. - z^^M 
Elik, l)a.(k, 1) + Elik, 2)a,(k, 2) = e^ - ^aJe^^M + ^ {S^ 
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Elik, l)a.(k, 1) + EUk, 2)a.(k, 2) = eL - ^.^^4.^ + ^ek. 
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With M =< S^ > and eko- = Ck — z^jjibB. Solving the above set of equations we obtain 
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This way of determining the spectral function and therefore the Green's function is known as the 
spectral density approach(SDA)^^. In this method, the only approximation is making an ansatz 



for S]ia{E). The rest of the calculation is exact. Therefore, the method is nonperturbative and 
is known to be more reliable in the strong coupling limit for the Hubbard model and Anderson 
model^^'^^. In materials that we are aiming to study also, it is the strong coupling limit since J 
is larger than the bandwidth. Having determined S'ko-(-E), the equation for the single-electron 
Green's function becomes 

G.,WH(c,„;eL)),^|:^^f|-i (12) 

Using the spectral theorem, the spin dependent occupation number is given by 

- ^ \^ Y" ( ^'^(^'■^') ^ n '^\ 

^'^ jY Z^ 2^ I /3(E,(k,i)-M) 111 ^^"^i 

k j=l \ / 

Here /3 = l/ksT. The chemical potential fi is fixed by the constraint 

n = ^^n(j = constant. (14) 

(T 

For the numerical evaluation of occupation numbers, the k-summation can be conveniently 
replaced by an integration over energy using a simple cubic (SC) density of states in a tight- 
binding approximation^^. In the next section we present the calculation of the static magnetic 
susceptibility. 

3 Static magnetic susceptibility 

In order to derive the magnetic phase diagram, one has to find out the magnetic ordering 
temperature Tc, which we find from the zeroes of the inverse paramagnetic susceptibility. The 
zero field static susceptibility of the itinerant electron subsystem is given by 



X(T) =^Bl]^a(^naj 



(15) 
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The local moment magnetization (S^) and itinerant electron polarization (n^ — ni) are mutually 
coupled and therefore they become critical for the same parameters, in particular, at the same 
temperature T. Therefore in the critical region, we can reasonably assume 
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The proportionality factor t] has to be still fixed. This we do later. A straightforward derivation 
of the itinerant-electron susceptibility eventually leads to the following expression: 

A{T) + B{T) \ 



xiT) = 2/i| (- 



-2r,J[B(T) + C{T)]j 
Where A{T), B(T) and C{T) are temperature dependent integrals which are given by 
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Where Q;(k, j) and i5(k,j) are the spectral weights and energies in the paramagnetic phase. 
The singularities are the solutions of the following equation: 



= 1 - 2r7J B{T) + C{T) 



T=Tc 



(22) 



The instabilities of the paramagnetic phase towards ferromagnetism are thus given by the 
solutions of the above equation. 



4 Results and discussion 

For numerical evaluation we use SC density of states (DOS)^° and the spin value •S' = |. We 
have chosen SC lattice density of states because it is one of the simplest with a peaked structure. 
In contrast, a rectangular or semi-elliptic DOS would be, as it is known for a long time, very 
unfavourable to magnetism. S* = | represents the situation for manganites. The width W 
of the Bloch band has been chosen to be leV. Our theory does not aim at a self-consistent 
determination of the local moment magnetization M =< S^ > but rather at the infiuence of 
the interband exchange on the band states. It is known from the D = oo limit of KLM, that the 
thermal behaviour of M is mean-field-like. We therefore consider M as a parameter for which 
we choose a Brillouin function. We have fixed /i from the requirement given in Eq.(14). In order 
to fix the factor rj which is introduced in Eq.(17). We study the temperature dependence of 
local magnetization (S^) / S and itinerant electron polarization (n^ — n^)/n in Fig.l. It is noted 
that both the magnetizations behave in a similar fashion. The similar behaviour of both the 
magnetizations particularly in the critical region where transition takes place from magnetic to 
non-magnetic phase suggests the following analogy: 
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Figure 1: Local magnetization (S^) / S (full line) and itinerant electron polarization (n^ 
n^)/{n^ + ^4.) (dotted line) with n = 0.4, S = 3/2 and J = 1. 



So that we assume rj = S/n. After fixing the value of /i for a given band occupancy n and 
using this value of rj, the paramagnetic susceptibility x of the itinerant electron system given 
in Eq.(18) is calculated as a function of temperature. 

In manganites the interesting physics is due to the dynamics of the electrons in the d-band. 
This band is split by the crystal field into two, namely the triply degenerate t2g-band and the 
doubly degenerate e^-band. In the undoped situation, the Mn is in the Mn^~^ and with four 
electrons in the d-band. The t2g band is occupied by three electrons with their spins aligned 
parallel due to strong Hund's rule coupling and this provides the localized spin S = 3/2. The 
Cg-band is occupied by one electron. Upon doping sufficiently Mn^+ state changes to Mn^^ 
state with three d-electrons. Then only t23-band is occupied and the Cg band is empty which 
means a hole is produced. Since in the present calculation n represents the number of itinerant 
electrons in the e^-band, the cases with n = 1 and represents the situation of Mn^^ state and 
Mn'^'^ state respectively. Therefore it is worthwhile to investigate the magnetic properties of 
the model system for given set of model parameters as a function of n. In order to investigate 
this, in Fig. 2 we have plotted the inverse susceptibility as a function of temperature for various 
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Figure 2: Temperature dependence of inverse susceptibility(dotted lines) and magnetization(full 
lines) for different values of band filling n. a) n = 0.05 b) n = 0.10 c) n = 0.40 and d) n = 0.44 
with S = 3/2 and J = 1. 



values of n. It is found that for various combinations of the model parameter values, the 
susceptibility follows Curie- Weiss behaviour. 

The magnetic ordering temperature Tc is determined from the zeros of inverse susceptibility. 
The n dependence of Tc is shown in Fig. 3. This figure clearly demonstrates that ferromagnetism 
does not exist if the band occupancy increases beyond a critical value and this critical value itself 
increases with the increase of the interband exchange J. Further, for any choice of the values 
of model parameters J and W, we did not find magnetism for n > 0.48. It is also noticed that 
the ferromagnetism is restricted to surprisingly low itinerant electron concentrations. It should 
be noted that the KLM used for the present investigation does not contain direct exchange 
between the localized moments. Therefore the collective ordering can only be mediated by 
conduction electron polarization. 

The steep increase of Tc for low band occupations, a rather distinct maximum (which shifts 
to the lower values of n with the increase of J) and then also a very rapid decrease to zero, qual- 



itatively mimics the experimental observation on diluted ferromagnetic semiconductors such as 
Gai^^MrixAs. A similar result is found with the "modified" RKKY theory^^'^^, interpolation 
self-energy ansatz method^^'^^ and also in another theoretical work^^ using the dynamical mean 
field approximation. Though we did not inspect, we speculate that the break down of the 
ferromagnetic phase is followed by some kind of antiferromagnetism. 
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Figure 3: Curie temperature Tc as a function of band filling n for various values of J and 
S = 3/2. 



Having determined the value of Tc, it is worthwhile to examine, particularly in the ferro- 
magnetic region, the quasiparticle density of states (QDOS) which we obtain from the single 
particle Greens function given in Eq.(12) through the relation. 
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The QDOS for various temperatures are displayed in Fig. 4. 

As is to be expected, the e^-band splits for each spin direction into two sub-bands centred 
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Figure 4: Quasiparticle density of states as a function of energy for various values of temper- 
ature. Full line for up-spin and dotted line for down-spin. The chemical potential is indicated 
by thin vertical line. Parameters: n = 0.40, S = 3/2 and J = 1. 

at ±JS. In the limit of large J the spectral weights and energies given in Eqns.(lO) and (11) 
reduce to 
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The separation of the bands and the spectral weights of these sub-bands however depend 
on T through the value of M. For example, at T = (M = S*), the spectral weight of the 
upper sub-band for f-states is zero. The reason for this is easy to understand. At T = 0, the 
local moment system is saturated. Therefore, for an '["-electron there is no chance to spin-flip 



by involving a corresponding spin-flip of the local moment system. That means, at T = 0, as 
far as the ^-electron is concerned, only the Ising part of Hgf operates resulting simply in a rigid 
shift of band without any splitting. The spectral weight of the J,-states in the lower sub-band is 
however not zero. This is because, for a 4,-electron, even at T = 0, spin-flip is possible. Further 
more, when a J,-electron flips its spin, it lands as an '["-electron. Therefore the QDOS of the 
J,-electron should be in the same energy region as that of f-electron. For a J,-electron there 
is another possibility. It can have repeated magnon emission and absorption. That is, in a 
sense, it propagates in the lattice dressed by a cloud of magnons. This is a stable quasiparticle, 
which we call the magnetic polaron. That is why for the 4,-spin there is a lower sub-band even 
at T = 0. However, at T = the width of the sub-band is infinitesimally small or the DOS 
is (5-function like and therefore is not shown in the flgure. Obviously, at T = there is no 
possibility of magnetic polaron for '["-electron. As T increases (Mdecreases from saturation), 
the spin flip processes are allowed for both spin directions and therefore the spectral weights in 
both sub-bands are non-zero for both spin directions, hi T = Tq (M = 0), the spectral weights 
of t- and 4,-states in the two sub-bands become equal as it should be. We note the asymmetry 
with respect to the center of the free band. This originates from the renormalization of the 
atomic levels by the H^f interaction^^. All these features can be understood from the QDOS 
for various values of the temperature. In plotting this we use the magnetization obtained from 
the Brillouin function for a given band occupancy. 

Further, it can be noted that at T = 0, for the lower sub-band, the width of the band, for 
J,-states is zero and for '["-states has twice the width of the free Bloch band and /i lies in the 
lower sub-band for n = 0.4. Therefore, the QDOS at fi has contribution only from '["-electrons. 
This situation in literature is described^^"^^ as a half-metal. The existence of half-metallic 
ferromagnets was flrst introduced by de Groot et al. In the case of doped manganites, half- 
metallicity was experimentally studied using low temperature resistivity data^'^ and by analysing 
the spin-resolved photo emission data^^. As the temperature is increased, the magnetization is 
reduced and the width of the band for states for '["-electrons decreases whereas for states with 
J,-electrons increases but the QDOS has contribution only from '["-states. When the temperature 
increases beyond certain critical value T^ (which we call the transition temperature from half- 
metal to metal), the width of the band for ^-states becomes sufficiently large so that the QDOS 
is contributed by 4,-states also. Thus one can see a transition from half-metal to metal as the 
temperature is varied. The critical temperature T^ naturally should and does depend on n. For 
further increase of temperature, beyond T^ the system undergoes transition from ferromagnetic 
to paramagnetic phase. It is found that Th depends on Tc and is less than Tc which is consistent 
with the experimental findings from spin resolved photo emission studies^ on LaQ^Sro^^MnOs. 
Since T^ and Tc both depend on n they can be varied by controlling the doping. In addition 



since DOS depend on M, they can be manipulated by changing the external magnetic field in 
place of temperature and thereby again T^ and T^ can be changed. Therefore one can drive the 
system from half-metal to metal by tuning the magnetic field or doping. Since the change in 
QDOS and the shift in ^ are continuous, the transition half-metal to metal is continuous. 

5 Conclusions 

KLM is used to study the possibility of collective magnetism in local-moment systems. SDA 
which is nonperturbative and is known to be reliable in the strong coupling limit for the Hub- 
bard model and the periodic Anderson model is used to solve KLM. Using this approach, 
single-particle Greens function and the occupation numbers are calculated for itinerant elec- 
trons in the e^-band. The zero field static susceptibility of the itinerant electrons which are 
exchange coupled to localized magnetic moment is derived. Then, by determining the magnetic 
ordering temperature from the singularities of the susceptibility, the magnetic phase diagram 
is constructed. It is found that the collective magnetism is possible only for low density of 
itinerant electrons. The QDOS of the e^-band is studied as a function of temperature which 
shows the possibility of half-metal to metal transition in the region where the magnetic order 
is present. The self-energy obtained in this approach is real which is an artifact of the ansatz 
made for the spectral function. Therefore the effect of quasiparticle damping is beyond the 
scope of the present investigation. 
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